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ABSTRACT 

We report two spectropolarimetric observations of SN 2005hk, which is a close copy of the "very pecu- 
liar" SN 2002cx, showing low peak luminosity, slow decline, high ionization near peak and an unusually low 
expansion velocity of only about 7,000 km s _1 . Further to the data presented by Chornock et al., (2006), at 
-4 days before maximum, we present data of this object taken on 9 November 2005 (near maximum) and 23 
November (+ two weeks) that show the continuum and most of the spectral lines to be polarized at levels of 
about 0.2-0.3%. At both epochs the data corresponds to the Spectropolarimetric Type Dl. The general low 
level of line polarization suggests that the line forming regions for most species observed in the spectrum have 
a similar shape to that of the photosphere, which deviates from a spherical symmetry by < 10%. In comparison 
with spectropolarimetry of Type la and Core-collapse SNe at similar epochs, we find that the properties of 
SN 2005hk are most similar to those of Type la SNe. In particular, we find the low levels of continuum and 
line polarization to indicate that the explosion mechanism is approximately spherical, with homogeneous ejecta 
(unlike the chemically segregated ejecta of CCSNe). We discuss the possibility that SN 2005hk was the result 
of the pure deflagration of a white dwarf and note the issues concerning this interpretation. 
Subject headings: supernovae - spectropolarimetry 



1. INTRODUCTION 

The nature of the complex explosion mechanisms thought 
to be behind Supernovae (SNe), of all types, is expected to 
leave a signature in the resulting geometry of the explosion. 
The study of the polarized light arising from SN explosions 
provides a unique and powerful method for directly measur- 
ing asymmetries in SNe. The power of spectropolarimetry 
is demonstrated by its application to the study of the shapes 
of SNe at large distances, providing an effective spatial 
resolution far superior to any current or curr ently planned 
imaging technique (see Wang & Wheeler 2008, for a review). 
In SN photospheres, the predominant form of opacity is due 
to electron scattering, which is also a polarizing process. In a 
spherical symmetry, light originating from different quadrants 
of the photosphere have opposite but equal polarization com- 
ponents, leading to a zero net observed polarization. In the 
presence of asymmetries, the cancellation of the polarization 
components arising from different portions of the ejecta is 
incomplete, leading to an observable polarization whose 
magnitude is related to the degree of the asymmetry. In 
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addition to the polarization of continuum light, the absorption 
due to line forming species exterior to the photosphere 
cast shadows on the photosphere. If the photosphere is 
not completely covered, a polarization signature will be 
associated with wavelength-dependent absorption features, 
with the degree of polarization related to the amount and 
orientation of the coverage of the photosphere by the line 
forming region. 

Previous polarimetric studies of SNe have shown that Type la 
SNe and Core-collapse Supernovae (CCSNe) have different 
polarization properties. Type la SNe are observed to have 
generally low polarizations, indicative of a close-to-spherical 
symmetry. The degree of polarization of these events de- 
creases with time, as deeper layers of ejecta are revealed, 
implying the ionization/ density structure becomes more 
spherically symmetric with depth. In contrast, CCSNe have 
been observed with a wide variety of polarizations, generally 
larger than those observed for Type la SNe. In addition, the 
degree of polarization is observed to increase with time for 
CCSNe, suggesting the ejecta become more asymmetric at 
increasing depth. 

In the classical observing modes of photometry and spec- 
troscopy, Type la SNe at certain phases can be confused upon 
casual inspection with Type Ic SNe at other phases due to the 
common display of Fe II lines. Special care must be given 
when trying to determine the type and the physics of "pecu- 
liar" Type la and "peculiar" Type Ic SNe. The application of 
spectropolarimetry can provide additional constraints, both in 
the abstract sense of the Stokes parameters and in the inferred 
geometries, that can help separate events belonging to the 
different distinct classifications. 

One particular class of SNe, which have been previously 
classi fied as Type la SNe (albeit "the most peculiar Type 
la"; |Li et al.||2003"l >, are the so-called 2002cx-like SNe. It 
has recently been suggested that members of this class may 
be representative of Ty pe la SN events th at result from a 
pure-deflagration (see Branch et al.| |2004] and references 
therein) or instead may be peculiar, low-luminosity Type Ic 
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SNe ( |Foley et al.|2009)|Valenti et al.|2009 1. 
Here we report two epochs of spectropolarimetry of the 
2002cx-like SN 2005hk, acquired at approximately max- 
imum light and two weeks later using the Very Large 
Telescope (VLT) of the European Southern Observatory 
(ESQ) SN 2005hk was discovered in UGC 272 by | Burket| 
[^Li] ( |2005] l. The location of SN 2005hk, relative to its host 
galaxy, is shown as Fig. [Tj Early spectroscopy was presented 



by|Jha et all (j20"06l>, I Chornock et aT] ( |2006| l, |Phillips et aL 
( |2007[ ) and lSahu et al.| ( |2008[ ). They dentify the early high 



ionization (e.g. Fe III) species typical of SN 2002cx-like 
SNe, as well as the bright Type la SN 1991T (see also 
|Serduke et al.||20"05 1, and the low velocities characteristic 
of S N 2002cx (~ 60 00 km s _1 ). The spectra of |jha et al. 
and [Chornock et al.l showed intermediate-mass elements 
with weake r features than re gular SNe la, but stronger than 
SN 199 IT. Chornock et al. suggested that Ni and Co may 
be present in the early spectra and hence in the outermost 
layers of the ejecta. The red portions of the spectra show 
appreciable features of O I and Mg II between 7500 and 
9000A, whereas SN 199 IT, a possible analog, showed a 
basically featureless spectrum in that phase and wavelength 
range. Both SN 2002cx and SN 2005hk showed similar fine 
scale structure suggesting that the structure was systematic 
and hence not the product of random clumping in the ejecta 
of either event. Phill ips et al.| ( |2007[ ) present extensive 
u'g'r'i'BVRIJHK photometry and optical spectroscopy of 
SN 2005hk, including the spectral data of |Jha et al.| and 



Chorn ock et al.| Importantly, they note that the fight curves 
of these two events are very similar, with neither showing the 
typical secondary maxima in R and I displayed by normal 
SNe la, but light curves that are also distinctly different than 
under-luminous SNe la such as SN 199 lbg. In addition, the 
bolometric light curve, including Swif t UV data, shows a 
rather slow light curve (Ami 5 (B) = 1.56; [ Phillips et al. |2007] l, 
uncharacteristic of under-luminous SN la and similar to 
norm al Type la SNe over a magnitude brighter at maximum. 
Chorn ock et al.| (|2006| presented spectropolarimetry of 
SN 2005hk, at 4 days prior to our first observation. They 
determined that the continuum was polarized ~ 0.4% in the 
red, with the data represented by a dominant axis. They also 
identified a weak modulation of the polarization of the Fe III 
line at A5 129. 

In Section [2] we present our spectropolarimetric observations 
and our adopted reduction procedure. The results from 
these observations are presented in ^3] and these results are 
discussed in In ^5]we present our conclusions. 



2. OBSERVATIONS 

Spectropolarimetry of SN 2005hk was acquired using 
the ESO-VLT on 2005 Nov 9 and 23, with the FORS1 
instrument in polarization PMOS mode. A log of these 
observations is presented in Table fT] The data were reduced 
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following the scheme outlined in Maund et al. ( 2007a| l. At 



each epoch, the observations consisted of four exposures 
with the half-wavelength retarder plate positioned at four 
position angles 0°, 45°, 22.5°, and 67.5°, to derive the two 
linear Stokes parameters Q and U 9 . Both observations were 
conducted with the 300V grism, without an order separation 
filter, providing a spectral resolution of 12. 5A as measured 



from arc line calibration data. Telluric lines were removed 
from flux spectra in the standard manner. Flux standard stars 
were observed, with the full polarimetry optics in place, to 
calibrate the flux levels of the SN data and remove telluric 
features. The telluric features are not expected to affect the 
polarimetry, except to decrease the level of signal-to-noise 
across the absorption, provided there is no significant change 
in the strength of absorptions between exposures with the 
retarder plate rotated by 45°. 

The instrumental stability between the observations at the 
two epochs was assessed by calculating the difference 
between the i nstrumental signature corrections e Q and e v 
(Maund 2008). At both epochs the difference was found to 
be + 0.1%, implying a limit on the instrumental stability 
and the measurement of the Stokes parameters to 0.1%. This 
is consistent with previous reported measurements of the 
instrumental signature correct ions and l imits on stability of 
the VLT FORS instrument ( |Maund et al.|2009] >. 
The degree of polarization was corrected for bi ases due to 
observ ational noise using the equations derived in Wang et aT] 
(1997), and, additionally, using a Monte Carlo sim ulation of 
the FORS 1 instruments in the manner of Patat & Romaniello 
( |2006[ ) and [Maund | ( |2008[ l. The wavelength scales were 
corrected for the heliocentric recessional velocity of the host 
galaxy of 3895 km s _1 . 10 

3. OBSERVATIONAL RESULTS 

3.1. The Distinguishing Spectral Features 

The spectroscopic properties of SN 2005 hk have already 
been described in some de tail in other works (Ch ornock et al.| 
2006| Phillips et al. 2007 ). Here the prominent spectroscopic 
features, which characterise SN 2005hk at the epochs at which 
spectropolarimetry was acquired, will be briefly reviewed to 
orient the reader. The flux spectra from our two spectropolari- 
metric observations are presented as Fig. [2] 
Our data from 2005 Novemb er 9 are inter leaved with the 
coverage presented by Phillips et al. (|2007 ) who give spec- 
tra from 7 and 14 November that they assign to -3 and +4 
days with respect to B maximum which they determine to be 
JD2,453,685.1. This means our data are from very near B 
maximum, perhaps a day earlier. We confirm that the spectral 
feature s are very sim ilar to those of SN 20 02cx at the same 
epoch ( Li et al.||2003] >. Furthermore, Chornock et al. (2006) 
acquired a flux spectrum of SN 2005hk in the course of ac- 
quiring spectropolari metry at N ovember 5.4, 4 days prior to 
our first observation. Chorno ck et al.| performed line identifi- 
cations using the SYNOW 11 code for their spectrum. They 
found that the strongest lines were from Fe III with other 
lines of intermediate mass elements being present, but rela- 
tively weaker. They obtained a better fit in the blue by adding 
lines of Ni II, Co II, Ti II and C III. They used an excita- 
tion temperature of 10,000 K for the low-ionization stage el- 
ements and 15,000 K for the higher ionization stages, Fe III 
and Si III. Our flux spectrum from the 2005 November 09 ob- 
servation, despite being acquired 4 days later, requ ires similar 
identifications to those made by Chorn ock et al.| In the red 
portion of the spectrum, we identify the Ca II IR triplet, al- 
though the synthetic spectra of Chornoc k et al.| suggested that 
at the epoch of their observations this was mostly du e to O I 
and Mg II. The feature at ~ 7600A is identified by Phillips 



We adopt the notation o fWang& Wheeler|{2008| , where the normalized 
Stokes parameters are given by Q = Q/I and U = U II. 



http://nedwww.ipac.caltech.edu/ 
1 ' http://www.nhn.ou.edu/ parrent/synow.html 



3 



et al. (2007) as Fe II, while Chornock et al. (2006) suggested 
it was a mixture of O I and Mg II. Li et al. ( 2003 i observed 
the appearance of O I A111A after maximum in spectra of 
SN 2002cx, which suggests that at days for SN 2005hk this 
feature could be at least partially due to O I. 
Our spectrum from 2005 November 23 corresponds to an 
epoch two weeks after B maximum. Our data are very s imilar 
to spectra of that date presented by Phillips et al. ( 2001) , but 
our data have higher signal to noise and broader wavelength 
coverage. The spectrum of Phillips et a l.|2007| from Novem- 
ber 25 (+15 days) has similar wavelength coverage to ours of 
November 23 and remains very similar to our flux spectrum. 
At this epoch, the underlying continuum is cooler with a rel- 
ative paucity of flux in the blue compared to the spectrum of 
9 November. The spectrum is characterized by the growth in 
amplitude of small scale features. S YNOW was used by us, in 
conjunction with previously published line identifications, to 
determine line identifications specific to the spectrum at this 
epoch. The degree of blending of multiple lines is lower at 
the second epoch, due to the relative decrease in the expan- 
sion velocity. The truncation of the red-side of the O I A111A 
emission line profile at 7980 A may be evidence for Mg II, 
but it is located close to strong telluric absorption features. 
The identification of Si II /16355 at +14 days is debatable, and 
is likely to b e due to an increase in Fe II AA6456, 6518 ( |Li| 
eTaLl[2003] ). |Branch et aL[ ( [2004] ) similarly found that Si II 
was no longer required to fit spectra of SN 2002cx at a similar 
epoch. The photospheric velocity was estimated, using lines 
of Fe 11(42), to be ~ -7000km s _1 . 

At the second epoch, there are also two distinct absorption 
features that are likely to be caused by the Ca II IR triplet. 
The Ca II IR triplet was identifie d at 12 days after maxi- 
mum in the spectrum of SN 2002cx (Branch et al.|2004|). The 



two features might be due to partially resolving the triplet in 
this low velocity photosphere (-5650 km s ), or there could 
be two kinematically separate co mponents such a s those ob- 
served in normal Type la SNe (Mazzali et al. 2005 ). If the lat- 
ter, the absorption minima correspond to velocities of -8400 
and -5500 km s _1 , respectively, taking the mean wavelength 
of the Ca II blend to be 8579A. The bluer of these two compo- 
nents is moving somewhat more rapidly than the typical range 
identifie d for photosphe ric features for Type la SNe at similar 
epochs (Mazzali et al. 2005). Comparison with SN 2002cx 
at 12 days (Br anch et al.||2004] ) leads to the identification of 
the strong absorption observed at 5730A as Na I intrinsic to 
the SN, corresponding to a velocity of -8150km s _1 (which is 

similar to the high velocity Ca II feature). 

At the first e poch, there is a feature at ~ 7200A. |Chornock| 
|et al.| ( |2006| ) argue it is unlikely to be due to C II due to ab- 
sence of other stronger lines that should occur in the wave- 



length range of these observations (e.g. /16578). Rather Sahu 



et al.|(|2008|) suggest s the 7200A may arise from Fe II ,17308 
Branch~et al.|(|2004 ) could not identify C I features in spec- 



tra of SN 2002cx. At the second epoch, however, there is a 
suggestion of an absorption (blended with an Fe II line giving 
rise to a flat topped emission line profile for that line) which 
is consistent with the expected position of C II /16578 with a 
velocity of -5500 km s (similar to the low velocity compo- 
nent of Ca II). 

3.2. The Polarization ofSN 2005hk 

Spectropolarimetry of SN 2005hk at the two epochs is pre- 
sented in Figs. [3]and|4] 



The correct measurement of the intrinsic polarization of SNe 
requires the subtraction of the polarization contribution of the 
interstellar medium, both in the Galaxy and in the host. In 
the first instance, a general limit can be placed on the degree 
of the Interstellar Polarization (ISP) by the amount of redden- 
ing measured towards SN 2005hk. For Galactic -type dust the 
maximum degree of the ISP i s related to the redden ing by the 
limit pisp < 9 x E(B - V) ( |Serkowski et al.|1975| ). The de- 
gree of Galactic foreground reddening in the directi on of SN 
2005hk has been reported to be E(B - V) = 0.022 ( |Schlegel| 
|et al.|[l998] ). The limit on ISP arising in the foreground is 
TO%. 

The reddening towards SN 2005hk can be estimated, if we 
assume that the narrow Na I D absorption component arises 
in the interstellar medium (ISM) of the host. The equivalent 
width of the Na I D absorption feature at the first epoch was 
measured to be W,\ = 0.3 ± 0.02. This value i s smaller 
than the value reported by Chornock et al. ( |2006| ), for ob- 
servations 4 days prior to our first observations, although the 
difference between t he two measured values is ~ 2<x. We 
note [Chornock et al. do not quote errors on their determined 
equivalent width. It is not clear if this difference is within 
the mea surement errors o r if it reflects evolution triggered by 
the SN ( |Patat et al.|2"007t. We follow |Chornock et al.|p006l ), 
who adopt the results of |Munari & Zwitter| ( |1997[ ) lor evalu- 
ating the relationship between Na ID absorption and the red- 
dening. T his relationship i s based on Galactic-type extinction 
laws (e.g. Cardelli et al. 1989 ). It is, therefore, implicitly as- 
sumed that the reddening component, measured from the nar- 
row Na I absorption component in the host galaxy, follows a 
similar law. As the Na I D component lines are not resolved 
in our data, we also adop t a range of flux rati os between Na I 
Dl and D2 of 1.1 to 2 ( |Chornock et al.|2006| ). This limits the 
reddening of the host to 0.05 < E(B - V) < 0.07. The total 
reddening, for Galactic type reddening and ISP law, yields a 
constraint on the tota l ISP of p ISP < 0.8%. 
Howe lTet al.| ([2001 ) present a technique for the determina- 
tion of the ISP for specific wavelength ranges which have 
null intrinsic polarization. This assumes that, for these spe- 
cific wavelength ranges, there is sufficient overlap between 
groups of lines such that there is complete in trinsic depolar- 
ization over these wavelengths. Howell et al. suggest that the 



wavelength region of 4800-5600A is intrinsically depolarized, 
where line blanketing opacity dominates over electron scatter- 
ing opacity. Any polarization, observed over this wavelength 
range, is expected to be due to the ISP alone. For SN 2005hk, 
at both epochs, the lack of significant overlap and blending 
of lin es leads to the c omplex behaviour of the Stokes param- 
eters. Chornock et al. ( 2006 ) find that the wavelength range 



of 4000-4200A satisfies the criteria for intrinsic depolariza- 
tion, namely sig nificant line blanketing. The observation of 
Chornock et al. at November 5.4 is sufficiently close to our 
observation at the first epoch that we observed similar features 
in the polarization spectrum and the Stokes parameter spec- 
tra. Importantly, we too determine that this wavelength range 
provides the best approximation to a line blanketed regime. 
The measured Stokes parameters across this range at the first 
epoch are Q = 0.04% ± 0.03 and U = -0.32% ± 0.03. The 
uncertainties on these quantities are the combined systematic 
and statistical uncertainties. 

The spectral evolution of SN 2005hk by the second epoch, 
with the increase in the strengths of individual Mg II and 
Fe III lines in the range of 4000-4200A, leads to the pres- 
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ence of intrinsic polarization (Q = 0.12% + 0.06 and U 

0.47% ± 0.060). Instead, we identify the region of 4250- 
4450A as being spectroscopically similar to a line blanketed 
regime with no unblended line structure. Over this wave- 
length range the average values of the Stokes parameters are 
measured as Q = 0.09% ± 0.03 and U = -0.30% + 0.03. The 
measured Stokes parameters, over the depolarized wavelength 
range, agree with those determined at the first epoch. 
We adopt values of the ISP Stokes parameters of Qjsp = 
0.07% + 0.03 and U ISP = -0.31% + 0.03, for both the Galac- 
tic and host ISP components, which ar e approximatel y con- 
sistent with the values measured by Chorn ock et al.| (who 
calculate separate Galactic and host components). The value 
of the ISP relative to the observed Stokes parameters is shown 
on Figs. [3] and [4] These values correspond to a polarization 
of pis p = 0.32% ± 0.03 with an angle of 9 ISP = 141° ± 3. 
The orientation of the ISP component is shown on Fig. [T 
As previously observed for other SNe (Mau nd et al.| [2007b~ 



Patat et al.||2009| |2010| i, the ISP component is aligned with 



the spiral arm pattern at the SN location, due to alignment of 
magnetic field lines and, hence, dust grains along spiral arms 
( |Scarrottetal.|1987| ). 

Given the limited wavelength range over which it can be de- 
termined, it is not clear if and how the ISP component varies 
with wavelength. Although a standard Serkowski-type law 
could be applied, it is unclear what the appropriate values of 
A max and p max would be. Furthermore, it has been observed 
for a number of SNe with significant ISPs that the ISP does 
not necessarily even follow a Galactic-Serko wski type law o f 
wavelength dependence (e.g. SN2006X; |Patat et al.||2009) . 
Given the small value of the ISP, we subtract, therefore, a 
constant ISP component from the observed data, to recover 
the intrinsic polarization of the SN (as shown by Pq in Figs. [3] 
andffl), under the assumption that any wavelength dependence 
of the ISP is smaller than the uncertainties derived above. 
On the Stokes Q — U plane (see Figs. Blandly, the ISP com- 
ponent is located close to the centroid of the SN 2005hk data. 
The calculated weighted centroids for the ISP-corrected data 

(g, TJ) are listed in Table [il 

After correction for the ISP, the intrinsic polarization of 
SN 2005hk is observed to be extremely low across the ob- 
served wavelength range. The degree of polarization of a 
number of features observed to be relatively strongly polar- 
ized in the uncorrected polarization spectrum, are no longer 
as polarized once the data have been corrected for the ISP. 
Utilising the data across the range 6600 - 7000A as repre- 
sentative of a region of pure continuum emission, the level 
of the continuum polarization is measured as 0.17 ± 0.05% 
and 0.19 ± 0.05% at the first and second epochs, respec- 
tively. At the second epoch the assumption that this wave- 
length range represents a true continuum region is most likely 
invalid, due to the presence of a number of minor, broad un- 
dulations which are most likely due to the lines (rather than 
incorrect flatfielding of the flux spectrum). This wavelength 
range avoids any polarization associated with the feature at 
7200 A, and is less than th e continuum polarization measured 
by |ChornocklM~aLl ( (2006) at 2005 Nov 5.4 of 0.36%. This 
level of polarization is con sistent with a s pheroidal configura- 
tion with axial ratio > 0.9 dHdflich|1991|). 



The significance of the polari zation signal is given by 
(S/N) P = P/o-p = y/Nj2P(S/N) ( |Patat & Romaniello|2006 1. 
At both epochs there are no features, at either 15 A or 30A 



binning, which are detected in the intrinsic polarization spec- 
trum Pq at > 3cr and only a sm all number of bins detected at 
> 2cr. |Chornock et al.| (2006) argue that in their data from 
5 November, only the line of Fe III /15129 is polarized and 
at a rather small level. On inspection, their data suggest that 
Fe III /14404 is also somewhat polarized. There is another fea- 
ture of similar amplitude showing a minimum at /14800 in the 
polarized flux that has no strong component in the total flux 
spectrum, but is likely to correspond to Si II A505 1 . 
Our data at the first epoch also show possible modulation in 
the polarization spectrum associated with Fe III /15129, as 
well as Si II /16355 and the 7200A feature. At redder wave- 
lengths there are a number of peaks which are associated with 
telluric features, suggesting that these were variable over the 
course of the observations. The apparent significance of these 
features is therefore overestimated, as the possible variability 
of telluric features was not considered in the reduction pro- 
cess and not accounted for in the final error analysis. 
By the second epoch, there are possible peaks in the intrinsic 
polarization associated with the absorptions of the Fe II mul- 
tiplet 42, the intrinsic Na I feature and the absorption of Fe II 
A6147. 

The absence of a significant signal in the polarization spec- 
trum does not necessarily imply that the intri nsic polar ization 
of s pectral features is also not significant. Wang et al. ( 1997} 
and Leonard et al.| ([2002 ) studied the change in polarization 
AP across particular spectral features, by comparing the mod- 
ulation in the Stokes Q and U parameters from the approxi- 
mate "continuum" level and at the wavelength of absorption 
minimum, where: 

AP = ^(Q(A co „ t ) ~ Q {Aline)? + (U (A cont ) - U (A lim )) 2 (1) 

AP is independent of the choice of the ISP component. The 
values of AP for particular spectral features, at both epochs, 
are presented in Table [2] 

Even when assessed by modulation of the Stokes parameters, 
the data at both epochs show few significant features. At the 
first epoch, Ca II H&K, Fe II and, importantly, Si II /16355 
show significant modulation, whereas at the second epoch the 
significant modulation is associated with the, now stronger 
Fe II lines. The interpretation of the data presented in Table 
requires some caution, however, as the selection of A co „, an 
X-line, by eye, can lead to a bias in favour of selecting those 
wavelengths that yield the maximum change in the Stokes 
parameter s. The prin cipal differences between the data pre- 
sented by Chornock et al. at -4 days, compared to our data at 
days, is due to much larger S/N achieved in thei r observa- 
tion. We estimate that at 7200A, at 50A binning, |Chornock| 
|et al.| achieved a final S /N ~ 1400, whereas our observation 
reached S /N * 390. 

At the first epoch, two sets of two contiguous bins (at 30A) 
are observed at both the 7200A feature and immediately red- 
ward at 73 15 A. An obvious line feature is associated with 
the former spectral line in the flux spectrum. The latter fea- 
ture corresponds to a telluric feature that, although it should 
be equally cancelled between the two beams at each retarder 
plate position, may appear polarized due to variability in the 
telluric line between the pair of observations used to derive 
each Stokes parameter. 

On the Stokes Q - U plane (see Figs. [5] and Rjl, the data 
at both epochs assume an approximate elliptical distribution, 
with no major structure deviating from the main concentra- 
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tion of points. At the first epoch, the data have no obvious 
preferred direction. By the second epoch the data are more 
extended across the Q-U plane, although again with no major 
deviations from the main concentration of data points. From 
Fig. [6] it is also evident that, by the second epoch, the degree 
of polarization has become more wavelength dependent, with 
data at bluer wavelengths being concentrated at lower values 
of Q and U than the data at the red extreme of the observed 
wavelength range. Our lower S /N data does not show the 
wavelength dependent elongation to the same degree as ob- 
served at -4 days, for which a dominant axis is clear to the 
eye. The behaviour of the data on the Stokes Q - U plane, at 



both epochs, is quantified and discussed in Section 3.3 



3.3. Decomposition of the Stokes parameters 

The determination of the dominant axis, and the de- 
composition of the Stokes vector, permits a correction for 
the random position angle at which target SNe have been 
observed, as well as determination of the principal axis 
of symmetry of the data, the relative orientations of the 
different line forming species and, ultimately, the underlying 
geometry. The dominant axis forms a line across the Stokes 
Q-U plane, and the decomposition process involves 
rotating the data so that the principal rotated Stokes axis is 
aligned with the axis of greatest variance for the data on 
the Stokes plane. The polarization signal is measured in the 
two orthogonal directions, which are effectively the "rotated 
Stokes parameters", along the dominant axis (Pd) and along 
the orthogonal axis (P ). A polarization component in the 
orthogonal direction permits the identification of deviations 
from the single axial symmetry represented by the dominant 
axis. 

Here we explore two complementary techniques for deter- 
mining the dominant axis: least-squares fitting of a straight 
line and a principal components analysis of the data on the 
Stokes Q-U plane. Importantly, these two different analyses 
should give identical measures of the dominant axis, but 
provide different parametrisation of the data on the Stokes 
Q-U plane. This latter point is particularly crucial for 
applying the SN spectro polarimetric classification scheme of 
Wang & Whee!er|p008l). 



3.3.1. Least-squares fit to the data 

Straight lines, of the form U = a + f3Q, were fit to the 
data of both epochs, using a weighted least-squares fitting 
technique, where the da ta were we ighted according to the 
errors in both Q and U ( |Press et al.|1992| ). In order to confirm 
the validity of this technique, the form of the straight line was 
reversed and the best fit was calculated for Q = 1//3U - a//3, 
giving identical results. The value of /3 yields the angle of 
the dominant axis with respect to the Stokes Q-axis, through 
which the data are to be rotated to decompose the Stokes 
vector onto the dominant and orthogonal axes. For a single 
axial symmetry, the dominant axis will pass through the 
origin of the Stokes Q-U plane. In instances where the data 
are offset from the origin (e.g. by a wavelength-independent 
continuum polarization component) the offset will be par- 
tially manifested in the measured a parameter. 
Dominant axes were calculated for the data at both epochs 
over two different wavelength ranges: 3700 to 8600A and 
4000-7000A. The former range wavelength encompasses the 
entire wavelength range of the observed data, whereas the 



latter excludes regimes with poor levels of signal-to-noise 
(S /N < 200). The regions of low signal-to-noise arise at the 
blue and red extremes of the data, due to line-blanketing and 
the blue spectral energy distribution significantly reducing 
the observed flux in the blue and red, respectively, and the 
declining detector response (in the red). 
The results of the least-squares analysis and the properties of 
the determined dominant axes are presented in Table [3] The 
values of a, ft, and the polarization angle of the dominant axis 
(Qdom - jarctan/J) were found to be relatively insensitive, 
within the uncertainties, to the choice of the wavelength 
range over which the dominant axis was calculated. The 
dominant axes are shown on the Stokes Q-U plane on Figs. 
[5] and [6] The large values of x 1 for the straight line fit, at 
both epochs, indicate that the data are poorly described by a 
straight line, with a significant polarization component in the 
orthogonal direction. At the first epoch, simple inspection of 
the data on the Stokes Q-U plane would suggest that the 
data are evenly distr ibuted in angle about the centroid (as is 
illustrated in Section 3.3.2 1. At the second epoch, a dominant 
axis is more clearly observed. 



3.3.2. Principal Components Analysis 

In order to better understand and quantify the distribution 
of the data on the Q-U Stokes plane and, in particular, the 
degree of the orthogonal polarization component, a principal 
component analysis of the data was conducted. 
The elements of the 2x2 weighted covariance matrix take the 
form: 



cov w (Xj, X k ) 



2; M~: Mi { X iJ ~ x j) { x «< ~ x >) N 



2/ 



_J 1_ 

AX, AX :l 



N - I 



(2) 



where /, k = 0, 1, Xq = Q, X\ — U (see Bevington & Robin- 
son |2003) , i = 1...N and the weighted average of the Stokes 

3" 



parameters are Xj - -^j > as shown in Table 



over the N wavelength bins. The eigenvectors and eigenval- 
ues of the covariance matrix provide an alternative measure 
of the angles of the dominant and orthogonal axes (9 e /2) and 
the ratio of the degrees of polarization in the dominant and 
orthogonal directions (an axial ratio for the resulting ellipse 
b/a), respectively. This latter quantity can be used to describe 
an ellipse on the Stokes Q-U plane which contains the data. 
The results of the principal component analysis are presented 
in Table [3] and corresponding ellipses are presented on Figs. 
[5] and [6] ^s expected, the polarization angle of the dominant 
axis determined using this analysis agrees, within the uncer- 
tainties, with the result from the least-squares fit. 
The non-zero values of b/a determined from the principal 
components analysis show that the data at the two epochs 
are inconsistent with both straight lines and null polarization. 
This suggests that, although the polarization measured at both 
epochs is low, we have measured a real polarization signal for 
SN 2005hk. The ratio of the dominant and orthogonal compo- 
nents measured at the second epoch is smaller than that mea- 
sured at the first, implying that at the first epoch the ejecta 
were approximately round and, with time, departed more sig- 
nificantly from the spherical symmetry. 
This interpretation requires some caution, however, as 



Chornock et al. observed data that corresponded to a dom- 



inant axis four days prior to our first observation, which we 
do not see in our data at the first epoch. This may reflect the 
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low S IN of our first observation, rather than an "oscillating" 
evolution in shape of the structure of the ejecta; with the ejecta 
appearing approximately spherical at maximum but departing 
from a spherical symmetry both before and after maximum 
light. 



3.3.3. The data on the dominant and orthogonal axes 

The data were rotated onto the dominant and orthogonal 
axes, using the angles determined from the least-squares fit 
analysis (for the limited wavelength range of 4000-7000A). 
The rotated data are presented as Figs. [7] and [8] Despite 
the apparent increase in the polarization angle from the first 
to second epochs, the determined polarization angles for the 
dominant axes are different fro m the angle measured at -4 
days, for which Chornock et al. determined O^om - 35.6°. 
At the first epoch, the dominant and orthogonal polarization 
components are similar to Stokes Q and U, due to the small 
angle through which the data were rotated. The polarization 
modulations associated with Ca II, Fe III, Si II are predomi- 
nantly observed in the dominant polarization component, and 
over the range 4000-7000A the deviations observed in the 
orthogonal polarization spectrum a re consistent w ith the ex- 
pected noise level. Similar to the |Chornock et al.| observa- 



tions at -4 days, an important exception is the 7200A ab- 
sorption feature (and possibly Fe III /15129) which has both 
domin ant and or thogonal polarization components (although 
|Chornock et al.| refer to the rotated Stokes parameters q^s p 
and ursp, these are equivalent to our dominant and orthog- 
onal axes). At our second epoch, the majority of features, 
associated with Fe II, Na I and, particularly, the 7200A fea- 
ture show some modulation in both the domi nant and o rthog- 
onal polarization components. While Chornock et al.| ob- 
served an obvious wavelength dependence in the dominant 
polarization component, with increasing Pj at redder wave- 
lengths (which we estimate to be 0.1% per 1000 A), 4 days 
later we found the wavelength dependence to be shallower 
0.05% ± 0.01 per 1000A. At +14 days we measure a similar 
wavelength dependence 0.07% + 0.01 per 1000A. This wave- 
length dependence is clearly seen on Fig. [6] The change may 
not be evidence for evolution in the wavelength dependence 
of the dominant polarizat ion component, in itself, as the ob- 
servations of Chornock et al. were acquired with higher S /N. 
Such wavelength dependence of the polarization, incr easing 
at redder w avelengths was observed for SN 1999by (Howell 
et al |2 001 ). At the temperatures determined for SN 2005hk 



and for Fe-rich matter, the opacity at the photosphere is dom- 
inated by b lends of line transitions for A < 6000A (see Fig. 
2 of [Hoflich et al.|1993 1. At redder wavelengths, the opacity 
is primarily due to electron scattering. The continuum po- 
larization may indicate an aspherical Thomson scattering re- 
gion, e.g. by a global asymm etry induced by th e rotation of 
the White Dwarf progenitor (Howel let al.|200f[ > or off-center 
delayed detonations (Ho flich et al.||2006[ ). Because multiple 
line scatterings are effectively depolarizing, the observation 
of increasing polarization with wavelength may be understood 
in terms of the decrease of the role of line opacity at redder 
wavelengths. If the dominant axis corresponds to a principal 
axial symmetry, then the line forming regions of the species 
observed in the flux spectrum follow the same axial symme- 
try as that of the continuum at the first epoch, but not at the 
second epoch. 



3.4. The Inferred Ejecta Geometry 

On the basis of the continuum polarizations measured at 
and 14 days, any asymmetry in the shape of the photosphere 
(assu ming a spheroi dal configuration) requires an axial ratio 
> 0.9 (Hoflich 1991). A similar limit on the asphericity of the 
photosphere is achieved if one uses the hig her value of con - 
tinuum polarization (0.36%) determined by |Chornock et al.] 
The general level of line polarization is low, suggesting 
that the line forming regions cover an approximately spher- 
ically symmetric photosphere homogeneously. Limits on the 
amount of polarization expected to be observed can be calcu- 
lated, given a level of continuum polarization and line depth, 
under the assumption that the line forming region blocks only 
unpolarized light from the photosphere, while not blocking 
polarized light arising from the limb (i.e. that the photo- 
sphere can be separated into discrete polarized and unpolar- 
ized zones): 



l conl 
'line 



(3) 



where I cont and are the observed fluxes in the continuum 
and in the absorption line minimum, respectively (i.e. such 
that the line depth is I cont - In ne - this convention is chosen as 
the observe d polarization ar i ses fro m the light that is not ab- 
sorbed). As Chornoc k et al.| (|2006) calculated for their obser- 
vation of SN 2005hk, we also calculate the theoretical limits 
on the line polarization for measured absorption line depths 
in our observations. 

At the first epoch, the maxim um polar ization for the Si n 
/16355 line is only 0.14%. As |Chornock et aL suggest, that 
the lines are not strongly polarized is commensurate with their 
apparent weakness. This is true at the first epoch, where line 
strengths of only > 0.7 are observed. At the second epoch, 
however, absorption line depths of ~ 0.5 are observed, and 
with the increase in depth of these lines it might be expected 
that, should the line forming regions be asymmetric, that there 
would also be a commensurate increase in the polarization ob- 
served across these lines (see Fig. [9j. 

Eqn. [3] assumes that all of the polarization is due to the con- 
tinuum, which implicitly assumes that the photosphere itself 
is asymmetric. In a number of cases, in particular for core- 
collapse SNe, t he classical limit for lin e polarization is ex- 
ceeded (see e.g. Ma und et al.|2 007a l 2009), for low values of 
continuum polarization which suggest approximately spheri- 
cal photospheres. Violations of this criterion can be explained 
by large scale structure, with a dissimilar geometry to that of 
the photosphere, or clumping, which may produce polariza- 
tion irrespective of the shape of the photosphere or the degree 
of continuum polarization. The degree of observed polariza- 
tion, for a clumpy line forming region, should depend on the 
size of the clumps: for a given line absorption depth, fewer 
larger clumps block larger areas with contiguous polarization 
vectors, whereas more smaller clumps can be distributed more 
evenly across the photosphere. 

In order to explore the possible role of clumping, and the 
expectation for the production of polarization across shallow 
lines, a Monte-Carlo toy model of a simple photosphere and 
obscuring clumped line forming region was constructed fol- 
lowing the scheme of Wang et al. (2007|. A circular contin- 
uum emitting photosphere was adopted, with the luminosity 
as a function of projected radius r p across a photosphere of 
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radius R p i, set to: 



/(r p )//(0) = l-k 



1 



(4) 



with limb darkening modelled using the coefficient k, which 
for k = yields no limb darkening. For a three-dimensional 
spherical photosphere, this would correspond to to a limb 
darkening law of I(r p = R ph sin 6)1 7(0) = (1 - k (1 - cos 6)). 
The radial dependence of the polarization was assumed to fol- 
low a quadratic form p cc (r p /R p h) 2 , with the constant of pro- 
portionality chosen to yield 15 % polarization for ligh t emitted 
at the limb of the photosphere ( Chandrasekhar 1960 ). Opaque 
circular clumps were assumed, with t >> 1, such that all 
photospheric flux from the region underlying the clump was 
blocked. For each iteration of the model, N clumps of size 
Rdump/Rph were randomly positioned across the photosphere 
and the corresponding flux of the uncovered regions, and their 
polarization components, were summed. Clumps were per- 
mitted to overlap and only partially cover the photosphere. 
For each size of clumps, 3000 simulations were conducted. 
This permits a comparison of the expected polarization for 
a given absorption line depth, assuming a particular clump 
size in the line forming region. The results of this model (for 
RdumplRph = 0.1, 0.2, 0.3 and 0.4) are shown as Fig. [9] 
For some of the models, zero polarization was achieved for 
uniform coverage of the photosphere by the clumps, however 
the difference between the average polarizations expected (as 
shown in Fig. |9]l and null polarization was > 1 cr for all the cal- 
culated models? These models suggest that, for a clumpy line 
forming region, the clump scale for SN 2005hk is very small 
(< Q.lRpf,) to produce lines of the depths observed, in the 
flux spectra at both epochs, with low polarization. The uni- 
formly low polarization observed for lines arising from differ- 
ent species suggest that the ejecta are generally mixed, with 
no separate structures on the plane of the sky for different el- 
ements or ionization stages(e.g Fe II and Fe III). The element 
responsible for the absorption at 7200A may be an excep- 
tion, with a line forming region that is more aspherical than 
other species. The data for SN 200 5hk, at the two epoc hs pre- 
sented here and those presented by Chomock et al. (2006), are 
consistent with an approximately spherical symmetric pho- 
tosphere surrounded by similar line forming regions for all 
species. 

The principal assumption of the model presented above is 
that the line forming r egion and the photosphere are sepa- 
rate. As discussed in 5 3.3.3 the opacity at the photosphere 
(where the quasi-continuum forms) is not due to electron scat- 
tering alone, but rather includes a significant wavelength de- 
pendent contribution from atomic transitions; this implies that 
line forming regions are not, in reality, separate entities except 
in the most high-velocity regions such as those discussed by 
|Kasenet al.| ( |2003] i. 

The size of the clumps and their relation to "large-scale" 
structure within the ejecta, leads to an interesting problem of 
semantics: since very large clumps (R c i„ m p/R p h = 0.4) may 
be considered large scale structure and very small clumps 
(Rciump/Rph - 0.1) tend towards constituting an effectively 
homogeneous medium. Large opacities require that a large 
fraction of the photosphere, the last Thomson scattering dis- 
tance, are reprocessed by lines, mostly due to iron-group ele- 
ments (which will result in a significantly altered spectrum). 

4. DISCUSSION 



4.1. SN 2005hk as a Type la SN 

The class of SNe similar to SN 2002c x currently has 15 
members, including 2005hk and 2008ha (|Foley et al.|[2009). 



As dis cussed in detail by |Li et aL ( |2003[ ) and Phillip s~et al.| 

005hk were of especially 



(2007), both SN 2002cx and SN 



low velocity at the photosphere near maximum light, ~ 6000 
km s _1 , had early relatively high ionization, and low peak 
luminosity. Neither showed the typical secondary maxima 
in R and I of normal SN la, but were also dissimilar to 
previously studied under-luminous SN la that do not show 



the secondary maxima (e.g. |Filippenko et al. 1992 1. There 
is also extreme diversity within the class, with SN 2008ha 
being « 3 magnitudes fainter than SN 2002cx (Foley et al. 
2009). The light curve decline is rather slow, unlike previous 
under-luminous SN la that serve to define the light curve 
brightness/decline relation prominently used in cosmology 
studies. 



Li et al. (2003) pointed out that no published model of a 
SN la corresponds to all the behavior of SN 2002 cx, but 
suggested that a pulsating delayed detonation ( |Khokhlo"v 
1991 ) model might be worth considering further. |Jha et aE 



(2006) suggested that the nebular spectral characteristics that 
they determined might be consistent with a pure deflagration 
model since the models in the literature have insufficient 
energy to correspond to normal SNe la and hence tend to 
have smaller expansion velocities and p roduce less 56 Ni and 
hence tend to be dim. Phil lips et al.| ( |2007| ) suggest that 
the low expansion velocities, low peak luminosity, and the 
absence of secondary maxima in the NIR light cu rves may be 
reproduced by a 3D pur e-deflagration model (|Gamezo et al. 



|2004||Ropke et al.|2006| that produces O.25M ~oF* r Ni. 
A key issue with the hypothesis that this class of explosions 
represents pure deflagrations of white dwarfs is that the pure 
deflagration models leave a substantial amount of unburned 
C and O in the outer layers. This C and O should show up 
in the early phases and while O I All! A is clearly present at 
the second epoch (and possibly the first), the case for C II 
/16578 (and possibly for a contribution to the 7200A feature) 
is less clear but it is not conclusively absent. In addition, 
|Kozma et al. ( 2005 ) suggest that the late-time spectra of pure 
deflagration models show no similarity to that of SN 2002cx 
either, in particular the absence of strong features due to C 
and O. Jha et al.| ( |2006[ ) claim a marginal detection of O I in 
late-time nebular spectra of SN 2002cx, with the suggestion 
that more unbumed C and O may be containe d inside the 
dense, still optically thick core at late times. Sahu et al. 
claim that at early times O I /17774 is suppressed by high 
temperatures and, at late times, nebular O I lines are also 
suppressed by high temperatures. The lack of observation 
of O 1 ,17774 for SN 2005hk at very early times (before our 
observations ) is uncertain, due to poor spectral coverage at 



early ti mes ([Phillips et al.||2007| ) and low S/N (Sah u et al. 
2008). |Sahu et al.| ( |2008| ) are able to approximately model 
the photospheric spectrum of SN 2005hk, based on the 



W7 deflagra tion explosi on model ( |Nomoto et al. 1984). 
Importantly, Sahu et al. find that, given the conditions of 
the SN, their spectral synthesis suggests large amounts of 
unburned O. Their model also assumes spherical symmetry, 
for which they achieve a reasonable approximation to the 
line profiles (which are somewhat s ensitive to any depa rtures 
from spherical symmetry; see e.g. |Tanaka et al.||2009 for a 
counter example). The features remain too broad, are for- 
bidden rather than permitted, and show prominent emission 
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of [O I] /16300, 6364. High densities, as evidenced by the 
low expansion velocities and P Cygni profiles observed at 
late nebular epochs, might suppress the [O I] feature, but 
it is predicted in the models and so its absence is a definite 
const raint on the nature of SN 2005hk. 
Chornock et al. ( 2006| l consider the possibility that these 



axis. The observations of SN 2005hk show it to have a 



events may be related to normal SNe la, but viewed along the 
"hole" left by having the binary companion nearly directly 



along the line of sight (Marietta et al. 2000; Kasen et al. 
2004). Such an aspect might increase the prominence of 
high-ionization species, weaken absorption lines, reduce the 
photospheric velocity and induce some polarization due to 
the non-spherical geometry of the ejecta. One issue with this 
model for SN 2002cx and related explosions is that it would 
result in higher luminosities than when viewed from other 
aspect angles. Even if the underlying explosion were a sub- 
luminous event such as SN 1991bg, if the view were "down 
the hole," the result would be brighter than SN 1991bg. Both 
SN 2002cx and SN 2005hk were sublumi nous by a bout 1 
magnitude in all bands (except H; Phil lips et al.|2007| l. It is 
certainly not clear that this is consistent with a "hole" model. 
IChornock~e t al. 2006 comment that their measurement of 
a small but appreciable continuum polarization and weak 
line features was inconsistent with the prediction of |Kasenj 
et al. ( 2004[ ) of a continuum polarization of ~ 0.1 percent 
and strong line polarization peaks for an explosion viewed 
slightly off the axis of the hole. As the SN evolves, the 
relative size of the hole and the photosphere will change 
(in particular if the hole closes, which |Kasen et al.| |2004| 
expect to be a slow process if it occurs), which should lead 
to a commensurate change in the polarization properties. 
Importantly, we see only moderate change in the polarization, 
despite the significant changes in the flux spectrum, between 
the epochs at which SN 2005hk was observed. In addition, 
it is not at all clear how such a model would account for the 



late-tim e permitted line spectrum delineated by |Jha et ak 
(2006). Our observations do suggest some line polarization. 



is demands some asymmetry in the ejecta, consistent with 
a single axial symmetry, but does not spell out quantitative 

agreement with the "hole" model. 

An alternative model might be a "plume" model ( |Plewa| 



et al. 2004 Livne et al. 2005) in which a single plume (or 
very few plumes) reach the surface, spreading fresh 56 Ni. 
If the plume burned little, the expansion might be slow and 
leave high densities at late times. Our spectropolarimetry 
suggests that the outer layers of the ejecta are approximately 
homogeneously mixed, although stratified, which suggests 
mixing in a large number of smaller plumes (with size 
< Q-2R p i, at the first epoch). 

This class of objects thus presents many conundra. They 
show high excitation early on but otherwise have none of the 
properties of SNe la like SN 199 IT that show early Fe III. 
Their light curves are dim but broad, violating the "Phillips 
relation". They have very slow expansion velocities and may 
remain partially optically thick at nearly 300 days after the 
explosion. They are mildly asymmetric in the continuum, 
in the Fe III early on and in Si II near maximum. The 
polarization of SN 2005hk is completely compatible with the 
observed range of polarization properties for Type la SNe 
(see Table |4). The level of polarization of SN 2005hk is 
not the lowest detected for a SN, with SN 1996X (a normal 



Type la SN; Wang et al. 1997| l showing polarization just 
above statistical significance and with no apparent dominant 



spectropolarimetric type of Dl (Wang & Wheeler 20081, 
similar to most Type la SNe, indicating that there is a 
principal axis of symmetry to the ejecta. Importantly, there 
are no apparent loop like features, even for the strongest lines 
(which is particularly pertinent for the observations at the 
second epoch) which would indicate a significant departure 
from this main axial symmetry. Normal Type la SNe, such as 



SN 2001el dWang et al.|2003a] l and SN 2004S ( |Chornock & 
Filippenko 2008), show some structure in their ejecta princi- 
pally in the Si II, as w ell as the high velocity component o f 
the Ca n IR3 (see also |Wang et al.|2003a[ |Patat et al.| 2009 ). 
Conversely, there are Type la SNe, such as SN 2004dt ( | Wang| 



et al. 2006), whose line polarization is principally in the 
dominant polarization component which is consistent with a 
single axial symmetry. In the case of SN 2004dt, however, 
the observed line polarization was ~ 1%. In summary, the 
similarity between the spectropolarimetry of SN 2005hk 
and Type la SNe indicates that their respective explosion 
mechanisms a r e both approximately spherical (see TableQ. 
Wang et al. (20071 identified a correlation between the 
degree of polarization of Si II, at -5 days, with the decline 
rate (Am 15 (B)) of the light curve. At some wavelengths 
SN 2005hk is an outlier from the Phillips brightness-decline 
relationship, although it is consistent with th e relationship 
in the / and H bands ( jPhillips "al|[2007] l. While the 
polarization measurements of SN 2005hk is offset fr om the 
measu red relation derived for normal Type la SNe by |Wang| 
|et al.| it is not as significant an outlier as the subluminous 
Typela SNe 1999by and 2005ke (see Fig. [10 1. 



4.2. SN2005hk as a Core-Collapse SN 

|Jha et al.j ( |2"006] l, |Valenti et al] ( |2009> and |Foley et al.| 
( 2009 ) consider the possibility that events such as SN 2002cx 
and 2005hk are not related to Type la SNe, but may in- 
stead be a type of CCSN. Based on the photometric and 
spectroscopic similarity between SN 2008ha and SN 2005hk 
with low-luminosity Type IIP and Type Ic SNe (2005cs and 
2007gr, respective ly) and the occurren ce of these events in 
late-type galaxies, Valen tTet al.| ((2009 ) suggest 2002cx-like 
SNe are low-luminosity Type Ic SNe, p otentially associated 
with Gamma Ray Bursts (GRBs). Foley et al. ([2009 ) present 
additional events in this class, which are claimed to have taken 
place in lenticular SO galaxies, which may argue against mas- 
sive stars as the universal progenitors for this type of event. 
Spectropolarimetry provides additional parameter space in 
which to examine the similarities between different types of 
SNe. 



Valenti 



et al.| ( |2009] l and |Foley et al.| ( |2009b note the spectro- 
umTIaritytetween SN 2005hk and SN 2008ha at +15 



scopic similarity 
and +8 days (after optical maximum), respectively, as well as 
their similarity with the subluminous Type IIP SN 2005cs at 



+44 days, w hilst it was still on the photometric plateau ( Pas- 
torello et al.|2009 1. A t such times Type II P SNe have gener- 



ally low polarizations ( |Leonard et al.|2006] l, seemingly similar 
to that observed by us for SN 2005hk. From a spectropolari- 
metric perspective, however, comparison of SN 2005hk with 
SN 2005cs whilst on the plateau is inappropriate. During the 
plateau phase, in Type IIP SNe, the photosphere is located 
in the hydrogen envelope. Due to the low ionization energy 
of H, the shape of the photosphere is approximately spher- 
ical. Once the photosphere has descended below the H en- 
velope, into the He layer (He being harder to ionize and ex- 
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cite), the shape of the photosphere more closely follows the 
shape of the und erlying exciting material. Despite an asym- 
metric explosion, Leonard et al. (2006) observed that during 
the plateau SN 2004dj had a low polarization, which increased 
significantly once the plateau phase had ended. Given the H- 
deficient nature of SN 2002cx-like SNe, a more appropriate 
comparison is with the H-deficient CCSNe, which have been 
observed (as listed in Table R} to have significant continuum 
and line polarization for different species indicating the role 
of highly asymmetric explosion mechanisms. Highly asym- 
metric interiors have also been observed for the Type IIP SNe 
2001dh and 2003gd (Maund et al., 2010, in prep.), which sug- 
gest that an asymmetric explosion mechanism is common to 
all these CCSN events. 

It is interesting to note another ongoing controversy involv- 
ing SNe la and SNe Ic. SN 2002ic was revealed by |Hamuy| 
et al. ( 2003 i to resemble a SN la embedded in a hydrogen- 
rich medium. |Hamuy et al. interpreted this as a SN la 



explod ing with an evolved binary companion. Wang et al. 
(2004) snowed that the hydrogen component was strongly 
polarized and argued that the geometry was consistent with 
an extended, clumpy, disk-like structure such as observed in 
proto-p lanetary nebula stars. Ham uy et al.| and | Wang et al.| 
( 2004 ) identified other events of this category that were orig- 
inally identified as Type Iln events, but which, on closer in- 
spection and comparison with SN 2002i c, revealed an under- 
lying Type I spectrum. More recently, Benetti et al. (2006) 
have challenged this interpretation, arguing that the underly- 
ing supernova in SN 2002ic more closely matched a SN Ic 
than a SN la. 

In summary, the nature of the polarization of SN 2005hk is in- 
consistent with the polarization measured for CCSNe, in par- 
ticular those of Type Ibc which have been generally observed 
to have very large polarization (especially line polarization) 
that increases with time. 

5. CONCLUSIONS 

SN 2005hk has been observed at two epochs, correspond- 
ing to maximum light and two weeks after maximum. Af- 
ter subtraction of the ISP component, the SN is observed to 
have a very low to null intrinsic polarization. A modulation 
in the Stokes parameters is observed across the Si II /16355 
line being the most significant feature, at the first epoch, with 
AP = 0.54 + 0.12%. 

At both epochs, a dominant axis can be fit to the data, al- 
though the data are poorly described by simple straight line 
fits on the Stokes Q - U plane. While there is modulation in 
the dominant polarization component, at the first epoch, that is 
associated with line features in the flux spectrum, the absence 
of similar modulations in the orthogonal polarization compo- 



nent implies that the line forming regions of the species in the 
flux spectrum share the same axial symmetry as that of the 
photosphere. The data at both epochs is classified as being of 
Spectropolarimetric Type Dl. 

In comparison with previously pu blished high S/ N observa- 
tions of SN 2005hk at -4 days ( |Chornock et al.||2006[ ), we 
observe a possible rotation of the polarization angle of the 
dominant axis, but it is unclear (due to the differences in S/N) 
if this represents a real rotation of the principal axis of sym- 
metry with depth into the ejecta. 

The low level of continuum polarization limits the asymme- 
try of the photosphere, at both epochs, to < 10%. The low 
line polarization is suggestive that the line forming regions are 
mixed, that different species occupy similar locations within 
the ejecta, and that they have a similar geometry to that of the 
photosphere. This is similar to the general behaviour of Type 
la SNe, and different from the chemically segregated ejecta 
observed for CCSNe. 

We conclude that the spectropolarimetry of SN 2005hk is not 
inconsistent with some of the expectations of the deflagration 
explosion model. We find that, in order to preserve the ap- 
proximate spherical symmetry of the photosphere and the line 
forming regions, any plume-like structures would be required 
to have a scale, in the plane of the sky, of < Q.2R p i, at both 
epochs. This would give the overall appearance of homoge- 
neously mixed ejecta containing both the products of nucle- 
osynthesis (Fe, Ni, Co) and unburned O at all depths. 
Spectral similarities between SN 2002cx-like events and low- 
luminosity Type IIP SNe are likely to be merely superfi- 
cial. The differences in the observed spectropolarimetry of 
SN 2005hk with CCSNe suggests, from a geometry perspec- 
tive, that these events have different underlying explosion 
mechanisms. Unfortunately, with spectropolarimetry of only 
one SN from the 2002cx-like class it is not possible to deter- 
mine if these results are applicable to the entire class, in par- 
ticular to the most peculiar member of this class, SN 2008ha. 
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TABLE 2 

Polarization changes across spectral features 



Species A cont (A) An„ e (A) AP 

2005 Nov 09 

Ca H H &K 3780 3960 

Fell 4225 4315 

Sill 6120 6270 

Fell? 7070 7215 

Ca II IR3 8340 8490 





2005 Nov 23 




Ca II H &K 


3750 


3930 


0.49 ± 0.23% 


Fe II 


4017 


4077 


0.54 ±0.19% 


Fell 


4610 


4910 


0.35 ±0.12% 


Fe II 


5910 


6030 


0.5 ±0.11% 


Fe II? 


7130 


7220 


0.34 ±0.14% 



0.64 ±0.15% 
0.34 ±0.11% 
0.54 ±0.12% 
0.45 ±0.17% 
0.31 ±0.40% 



TABLE 3 
Best-fit Dominant Axes 



Least Squares Fit Analysis Principal Components Centroids 



Date 


a 


P 


X 2 /dof 




Axial Ratio 
b/a 


e /2 


Q 






U 


09 Nov 2005 A 


0.138 ±0.008 


0.077 ± 0.074 


347.2/154 


2.2 


t 7.7 


0.87 


3.7 


-0.037 ± 


0.016 


0.135 


±0.016 


09 Nov 2005 B 


0.138 ±0.009 


0.226 ± 0.084 


202.3/98 


6.4 


±7.9 


0.78 


8.7 


-0.049 ± 


0.016 


0.127 


±0.016 


23 Nov 2005 A 


0.075 ± 0.008 


0.545 ±0.051 


288.8/154 


14.3 


±5.0 


0.47 


14.0 


0.011 ± 


0.016 


0.081 


±0.015 


23 Nov 2005 s 


0.079 ± 0.010 


0.771 ±0.077 


160.9/98 


18.8 


±5.0 


0.42 


18.55 


-0.022 + 


0.016 


0.062 


±0.016 


A for data over the wavelength range 3700-8600A 



for data over the wavelength range 4000-7000A 



TABLE 4 
Comparison SN Spectropolarimetry 



SN 


Spectral 


Epoch ' 


SP 


Continuum 


Line 


Ref. 




Type 




Type" 1 " 


Pol. (%) 


Pol. (%) 




1996X 


la (normal) 


-7 


Nl 


< 0.3 


Si 11:0.6 ±0.15 


1 


1999by 


la (sub) 


-1 


Dl 


0.8 ±0.1 




2 


2001el 


la (normal) 


+2 


Dl+L(Ca) 


0.3 


Si:0.3;Ca:0.7 


3 


2003du 


la (normal) 


+ 18 


Dl 




0.2 


4 


2004S 


la (normal) 


+9 


D 


0.4% 


< 0.5 


5 


2004dt 


la (normal) 


-7 


Dl(Si+Mg) 




Si 11:2 


6 


2006X 


la (normal) 





Dl+L 


< 0.2 


Si 11:0.8 


7 


2005hk 


2002cx-like 


-4 


Dl 


0.4 




8 









Dl 


~ 0.2 




this work 






14 


Dl 


~ 0.2 




this work 


2002ap 


Ic (broad-lined) 





DL 


~ 0.3% 


Oi:4 


9 


2005bf 


lb (pec) 


-6* 


Dl,L(He+Ca) 


< 0.45% 


He I:1.3,Ca 11:4 


10 


2007gr 


Ic 


+21 


Dl+L(Ca) 


0.5% 


Ca 11:2.5 


11 


2008D 


lb 





Dl+L(He+Ca) 


< 0.4% 


He I:0.6,Ca 11:1.8 


12 






14 


Dl+L(He+Ca) 


< 0.4% 


HeI:l.l,CaII:2.5 


12 



' Relative to optical maximum; * SN Spectropolarimetric types see 
Wang & Wheeler 1 2008 1; * relative to t he sec ond maximum; References (1) 
Wang e t al.NlW7j; ( 2) lHowell et al.|<200"r); (3)|Wanget al.| 12003a); (4) 
Leonard ~ al.H2005>; (5) |(Jhornock_& Filippenk"o| |2U(J8( ; (6) |wang et al 
i'M6\; ( 7)|Patat et al.|j2CK)^t; ( 8 )|Chornock et aTl|2Tjn5) ; ( 9) | Wang et al. 

2003b I, Kawabata et al. 1 2002 1. Leonard et al. (20021; (10) Maund et al. 

2007aj; (1 l) |Tanaka et ai1 p308 1; ( 1 2) Maund et al. 1 2009 1. 
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Fig. 1. — SDSS r'-band image of the site of SN 2005hk, indicated by the 
filled circle, relative to its host galaxy UGC 272. The orientation of the de- 
termined ISP component, at PA = 141 ± 3°, is shown by the arrow. 




Fig. 2.— Flux spectra of SN 2005hk at 09 Nov 2005 and 23 Nov 2005, or 
approximately at B light curve maximum an d 14 days after maximum. Line 
identifications are based on those provided by Chornock et al. 1 2006 l, |PhiUips| 
|et al.|j2007} , |Branch et al.|(20()4 and fits to the spectra using the SYNUW 
code. The wavelength scale has been corrected for the recessional velocity of 
the host galaxy. 
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4000 5000 6000 7000 8000 



Wavelength (A) 

Fig. 3. — Spectropolarimetry of SN 2005hk on 2005 Nov. 9. The six pan- 
els (from top to bottom) give: the observed flux spectrum, in logarithmic 
flux units (ergss -1 cirT 2 A~'); the observed polarization spectrum (P); the 
observed polarization spectrum corrected for the ISP component (Pq); the 
Stokes Q and U spectra uncorrected for the ISP (the value of the ISP is shown 
for each parameter by the grey line); and the polarization angle 8. Line iden- 
tifi cations are provided in t he top panel, and are b ased on the identifica tions 
by |Chornock et al.|p006) , |Branch et al.| fzOQ4l, Philli ps et aT]j2007) and 
comparison with synthetic spectra produced using the SYJNUW code. The 
data have been rebinned to 30A for clarity. The wavelength range selected 
for the determination of the ISP are shown by the heavy lines. 




Fig. 4.— 
Nov. 23. 



Wavelength (A) 

Same as Fig. [3] but for spectropolarimetry of SN 2005hk on 2005 
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Fig. 5.— Spectropolarimetry of SN 2005hk on 2005 Nov. 9, on the Stokes 
Q-U plane. The data have been binned to 30A for clarity, and have been cor- 
rected for the ISP component, which is indicated on the Stokes plane by the 
filled black circle, with the radius of the circle corresponding the measured 
uncertainties. The data are colour coded according to wavelength following 
the scheme of the colour bar. The origin of the Stokes plane is indicated by 
the cross (+). The black dashed circle represents the null polarization enve- 
lope, at 6000A. The heavy black line indicates the best fit to the "dominant 
axis." The ellipse shows the apparent axial ratio and orientation of the data, 
as determined using the Principal Components Analysis, with the size of the 
ellipse arbitrarily scaled. 



23 Nov 2005 




4000 5000 6000 7000 8000 

Wavelength (A) 



-1 -0.5 o 0.5 l 

Q(%) 



Fig. 6. — Same as Fig. |5]but for spectropolarimetry of SN 2005hk at 2005 
Nov. 23. 
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Fig. 7. — Spectropolarimetry of SN 2005hk on 2005 Nov. 9, projected onto 
the dominant axis (Pj) and the orthogonal axis (P ). The data have been 
binned to 30A for clarity, and have been corrected for the ISP component. 
The weighted average of the ISP, rotated onto the dominant and orthogonal 
axes, is indicated by the heavy grey lines. 
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Fig. 8. — Same as Fig. [7] but for spectropolarimetry of SN 2005hk on 2005 
Nov. 23. 
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Fig. 9. — Results of Monte-Carlo simulations of random obscuring spherical 
clumps in front of a spherical photosphere. The models are shown for clump 
scales ofRciwnp/Rph = 0.1 (blue □), 0.2 (red x), 0.3 (green o) and 0.4 (black 
*). Solid lines are for simulations with k = 0.5, while dashed lines indicate 
simulations with k = 0. The shaded area indicates the permitted values of 
line polarizati on, according to E qn. |5] for a continuum polarization of 0.36% 
determined by Chornock et al. (2005). Lower values of the continuum polar- 
ization (0.17% determined here) place even stricter upper limits on the scale 
of clumping. Minor fluctuations in the curves are due to the statistical na- 
ture of the Monte Carlo simulations; except for the case of Rclump/^ph = 0.4 
for Iiine/Icont > 0-8, where the results are dominated by the finite size of the 
"pixels" that compose the model photosphere. 




Fig. 10. — Degree of polarization across Si II /16355 as a function of light 
curve decline rate. Our measurements of SN 2005hk are shown in red, with 
the value of AP(Si II) indicated by the point with the heavy outline, data in 
black are ta ken from|Wang et al.H2007t, the measurement for SN 2006X 
(A) is from |Patat et al.| (2009), the meas uremen t for the luminous event 
SN 2009dc (□) is taken from Tanaka et al. (2010), and the measurement for 
SN 2()05ke (x) is taken from Patatetal., (2010, hi prep.). 



